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The Helmholtz free energy, A, of a rigid body is a 
function of temperature, and of the six homogeneous strain 
components. If the crystal is to be rigid, three inequalities 
must be satisfied for the derivatives of A with respect to 
the six strain components, for a regular (cubic) lattice. 
This enables one to limit the pressure-temperature range 
for which the crystal is stable. The violation of the condi- 
tion cy >0, that the crystal resist shearing, is interpreted 
as leading to melting. From a knowledge of the forces 


HERE exist many attempts to derive the- 

oretically the laws of melting. I mention 
Lindemann’s! formula, Griineisen’s? general ther- 
modynamics of isotropic solids, the theories of 
Braunbeck,? Raschevsky* and of Herzfeld and 
Goeppert-Mayer.’ The newest contribution to 
this problem is an ingenious paper of Lennard- 
Jones and Devonshire® who apply the “method 
of order and disorder’ invented by Bragg and 
Williams.’ 

All these papers use different. ‘‘criteria’’ of 
melting: e.g. Lindemann, the collision of a 
vibrating molecule with its neighbor; Raschev- 
sky, the instability of a particle under the action 
of its neighbors; Braunbeck, the instability of 
the vibration of two rigid simple lattices against 
one another; Herzfeld and Goeppert-Mayer, 
the minimum of pressure with respect to changes 
of volume; Lennard-Jones, the transition of the 
structure from order into disorder. 

In actual fact there can be no ambiguity in the 
definition of, or the criterion for, melting. The 
difference between a solid and a liquid is that the 
solid has elastic resistance against shearing stress 
while the liquid has not. Therefore, a theory of 
melting should consist of an investigation of the 
stability of a lattice under shearing stress. As far 
as I know, the only author who has made an 
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4N. v. Raschevsky, Zeits. f. Physik 40, 214 (1927). 
5K. Herzfeld and M. Goeppert-Mayer, Phys. Rev. 46, 
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between the molecules the phase integral, and therefore 
the free energy, may be calculated as a function of T, V, 
and the six strain components. The numerical calculations 
are carried out for a body-centered cubic lattice. The 
product of all the frequencies is calculated directly, so that 
the assumption that the Debye equation for the frequency 
distribution holds, is not necessary. The melting curve, 
pressure against temperature, is then determined. 


attempt to apply this natural definition of melt- 
ing is Brillouin.’ The short published account of 
his considerations indicates that he uses Debye’s 
approximation for calculating the frequencies 
of the crystal; this method is objectionable, not 
only because it neglects the short waves, but 
because it does not lead to a reduction of the 
observable quantities to atomic forces. 

Independently of this work of Brilloain, I 
have developed a method for treating thermo- 
dynamics of a crystal lattice in such a way that 
the formulae are valid over a wide range of tem- 
perature, in the hope that this range might in- 
clude the melting point. 

This theory is in some way the counterpart of 
Mayer’s rigorous treatment? of condensing gases 
where the liquid state was reached from the 
other side, at least in so far as the existence ofa 
condensation point could be proved with the 
methods of kinetic theory of gases. 

The chief point in our method is the clear dis- 
tinction between molecular variables and molar 
parameters as defined in statistical mechanics. 
The molecular variables are the generalized 
coordinates and conjugate momenta 4, 1; 
gx, po; +++ which are distributed at random ac- 
cording to statistical laws. The molar parameters 
@1, @2, --- describe the (macroscopic) external 
influences on the system (e.g. the volume of a 
gas). The energy of the system depends on both 

8L. Brillouin, Phys. Rev. 54, 916 (1938). 

9J. Mayer, in collaboration with Ph. G. Ackermann 
and S. F. Harrison, J. Chem. Phys. 15, 631 bare 
(1937); III, 6, 87; 1V, 6, 101 (1938). 
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391 (1938); B. Kahn and G. E. Uhlenbeck, Physica V 


4, 399 (1938). 
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kinds of variables, €(91, £1; G2, P23 °° * 341, G2, ** -). 
If it is given, we can construct the partition func- 
tion which for high temperatures is given by the 
law of classical statistical mechanics: 


QO(a1, G2, ***; r= |: ; » { dgsdprdae of: 


Kea, Pll q2*+33 ay+++)/kT. (1) 
Then the free energy of the system is 
A=—hkT log Q, (2) 


and all other thermodynamic properties are 
found by differentiation : 


the entropy S=-—0A/dT, 


3 
the energy E=A+TS, (3) 


the generalized forces corresponding to the 
molar parameters: 


F,= —0A/da,. (4) 


If we want to apply this prescription to a 
crystal lattice we must clearly state which are 
the molar parameters. For the sake of simplicity 
let us consider only a cubic lattice of the Bravais 
type (simple, face-centered or body-centered). 
As thermodynamics has only to do with states of 
equilibrium (from the macroscopic, or molar, 
standpoint) we can restrict our considerations to 
homogeneous deformations. Then the cubic cell 
(with the side a) becomes a general paral- 
lelepiped, described by three vectors A}, Ao, a3 
with six essential parameters (invariant with 
respect to rigid motions), viz. three lengths and 
three angles, or the six scalar products of the 
vectors with one another 


9 
@y°, Ao”, as’, @e-A3, Asi, a1ao. 


The shape of the cell, i.e., these six parameters, 
are given by macroscopic (homogeneous) de- 
formation and play the part of molar parameters. 
The actual positions of the particles can be 
described by displacement vectors relative to the 
deformed cell, u’, where the index / stands for 
three integers 1,, lo, 1, defining a lattice point. 
The energy is a function of all the u’ and the 
ay’, -+-@o-as, ---. If we develop it with respect 
to powers of the u! the linear terms vanish in 
consequence of the symmetry of the lattice, and 
the quadratic terms allow the introduction of 
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normal coordinates which can be used, instead 
of the u’, also in the higher terms. Then the 
integration in (1) can be performed, at least in 
principle, and leads to 


(5) 


It is evident that A is proportional to the num- 
ber N of the particles (if surface effects are 
neglected). We introduce, instead of the a;’, 
+++ @o-@3, +: the side of the cube of the original 
(undeformed) cell and the six dimensionless 
strain components* 


A=A(T, ay”, as, a3", @2°@3, @3°A1, Ar" a2). 


2x2= 11 = (a1? —a?) /a?, Y2= €23=€2°a3/a’, 


2y y= C22 = (as? —a’)/a’, Z,=€31=43'a;/a*, (6) 
22,=€33= (a3’—a’) /a’, Xy= €19= 81° a2/a?, 
and write 
A =Nf(T, a, Xz, Vy Zz, Vey Zz, a) (7) 


The function f depends on T and 7 geometrical 
variables, but only six of the latter are essential ; 
there exists a linear identity between the seven 
first derivatives of f, 


adf/da=(2xz#1)0f/dxz+ +--+ +2y,0f/dve+ +-:. 


By repeated differentiation one gets identi- 

ties between the higher derivatives. We write 

down the two equations which hold for the 

first and second derivatives in the case x,=0, 
- y,=0, °«*, 











ay Df aie Os 
a— =——_+—-+—, 

0a OX, OVy O28; 
anh of je o7f of 
a*— — a— = —-++—_+—- 
da da Ox,” Oy,? 82,2 

é°f a°f a°f 
+2 +2 +2 ~ Sy 
OY ,02, 02.0%, OX20Vy 


If we develop f for a given value of a into a 
power series of the xz, +--+ yz, ++: the result must 
have, for a cubic lattice, the form 


f=fothi:(xz+yy+22) 


+3 (fi 5 (x2 +y/+2,7) 
+ 2fia* (YySetSe2+x2Vy) 


+fas: (v2 +22+%,7)}+ ns (9) 
* We use two different notations, that of Voigt: xz, «++, 
Ye, ***, and that of Brillouin: e::, «++, é23, -+-—. It is 


essential to define the strain components for finite deforma- 


tions, as we have done here, not in the usual way for small 
displacements. 
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where fo, fi, fis, fiz, fas are functions of a and T; 
the relations (8) become 


adfo/da = 3fi, a0°fo/da® = 3} fitfiut2pfis} ; (10) 
the latter can also be written 
adf;/da= 2fi +fiit2fic. (10a) 


The density of free energy is A/V where V is 
the volume containing N particles. Let y be the 
number of cells a* corresponding to one particle, 
then V=yNa*. The value of y differs for the 
different cubic lattices; if a@ is always the pro- 
jection of the smallest distance on the side of 
the cube, we have for the simple cubic lattice: 
y=1; body-centered : y =4; face-centered: y =2. 
Now we get for the density of free energy 

ee. fo 


a *(a,+ +2.) 
V ya? ya’ 2 be 

+3 {€11(x2?+y,?+2,7) 

+2610(VySZe+23XH2tX2Vy) 


+ Caa(Ve?+227+%,7) } a eh ’ 


—p=(1/ya*)fi(a, T), * 
Cu=(1/ya*) fila, T), 
C12=(1/ya*) fie(a, T), 
C44= (1/ya*) fas(a, T). 


The derivatives of A/V with respect to xz, 
+++ y,, +++ are the negative components of the 
stress tensor: 


—X,= —pteixetCre(yy +2.) + “sa. 


(11) 
(12) 


where 


(13) 


: (14) 
= V,=Cayz+ Sie 


We see that the stress does not vanish for x,=0, 

-y,=0,--- but goes over into a uniform 
pressure p=X,=Y,=Z,. This must be so as 
the constant a of the cubic lattice, to which the 
stress X,,--- Y,,--: is applied, is not assumed 
to be that of equilibrium under zero pressure. 
Our theory corresponds exactly to the experi- 
mental conditions under which Bridgman and 
others!” have recently worked, namely a combi- 


10 P, W. Bridgman, Am. J. Sc. 36, 81 (1938). J. App. 


. 9, 517 (1938). ; . 
pha ee p. 07 (1939). Here a list of the literature 1s 
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nation of hydrostatic pressure and other stresses. 

The formula (12) is evidently the equation of 
state for the cubic crystal giving a as a function of 
pb, T, hence the thermic expansion for any 
pressure. 

C11, Ci2, Cag are the elasticity constants in 
Voigt’s notation for a given temperature and 
external hydrostatic pressure. 

The stability conditions, which guarantee that 
the quadratic terms in (11) are positive definite,* 
are 


3 
—=C)+2¢12.>0, 


K 


C11—Cy2 > 0,7 C4a > 0, (15) 


where x is the cubic compressibility. The identi- 
ties (10), (10a) allow us to write the equation of 
state and the expression for the compressibility 
in the form 


p= <a) — (a— — 2f:) = : 
3ya? da sok «=Oya®X\: Oa da 








(16) 


which can be used for checking calculations. 

We now discuss the different possibilities for 
instability of the lattice if temperature increases. 

I. The normal case will be that where the 
stability against shearing stress vanishes first, 
C44=0. According to our previous considerations 
this means melting. 

The melting curve is given by the two equations 


fila, T)=—ayp, — fas(a, T)=0. (17) 


The second one gives a(7), and if we introduce 
this in the first equation we get p(T). 

Il. It might happen that ¢11=Cc1 is satisfied 
first so that the two equations 


fila, T) a —a*yp, fila, ie) —fir(a, T) =0 (18) 


are satisfied. This means the transition of the 
solid into a state in which the stress is a hydro- 


* Our minimum problem is not quite of the usual type 
because the variables a, xz, -*:, Yz, *** are not inde- 
pendent ; x: +9,+2: is the relative change of volume which 
can also be expressed as 36a/a. The orthodox way of 
treating the problem is this: We consider only 6 variables, 
say d, Xz, Vy) Ye, Zz, Xy by eliminating z, with the help of 
Xetyyt2z-=0. Then %z, yy are independent just as 
Ve) Zz, Xy- The linear terms in (11) vanish, and the quadratic 
terms can be written in the form 


4(¢11— C12) ° {3(x2+yy)?+(x2—Yy)*} + 3cu(yP+222+x,?) ; 
this gives the second and third of the conditions (15). 
The first follows directly by developing with respect to a 
and using (16). 
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static pressure (as in a liquid) but nevertheless 
an elastic resistance against shearing survives. 
I suggest calling this state a “gel.” 

Dr. K. Fuchs has directed my attention to the 
fact that there are indications of the existence of 
such a state even for metals. Measurements of 
the elastic constants of the alkali metals Na and 
K made by Bender" have shown that the 
difference ¢11—¢12 is unusually small. Bidwell” 
has measured electric resistance and thermo- 
electric power of these metals at various temper- 
atures and found a change of slope at about 100° 
below the melting point indicating a transfor- 
mation; but a study of the structure with the 
help of x-rays showed no change of the lattice 
type, but only an “‘obliteration”’ of the structure. 

III. If c1:-+2¢12=0 is satisfied first so that 


tila, Tp) = —a* yp, fila, T)+2fi2(a, T)=0, (19) 


the lattice is completely unstable and has no 
cohesion at all. . 

This can be interpreted as sublimation; but a 
complete coincidence with the observed sub- 
limation curve can hardly be expected because 
the representation of the phase space by the 
normal coordinates of a lattice excludes just 
that part of it which is characteristic for the 
gaseous state, namely the domain corresponding 
to free particles. The curve (19) will give a 
finite value of T; for p—0 (just as the melting 
curve (17)) whereas the real sublimation curve 
runs into the point 7=0, p=0. Nevertheless I 
think that for high temperatures, near the 
melting curve, the formulae (19) might give a 
fair approximation of the sublimation curve. 

If these curves intersect there exists a triple 
point, the solution of the three equations: 


4Q. Bender, Ann. d. Physik 5, 359 (1939), gives the 
following experimental data for Na and K (in units of 10" 
dynes/cm?), to which are added the theoretical values for Li 
(K. Fuchs, Proc. Roy. Soc. A157, 444 (1936)). 





Cll —Ci2 C44 
Na 0.166 0.618 
K 0.087 0.263 
Li 0.341 1.329 


This gives for (Ci1— C12) /C4g the values 0.27, 0.33, 0.26 in 

tie be gt arbi In general the value of the ratio 

is about ten times larger, e.g. 2.7 f 2 

ets 3.5 for fluorite. ok os au 
2 C. C. Bidwell, Phys. Rev. 2 : 

(1926). ys. Rev. 23, 357 (1924); 27, 381 
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fila, T)= —a* yp, fas(a, T) =0, 
fula, T)+2fi2(a, T) =0. 
But it also is possible that these equations have 
no solution (for positive values of a, p, T); such 
a crystal would never melt but only sublime 
(as diamond, for all attainable pressures). 
Two other triple points can exist, with respect 


to the “gel” state. 
Entropy S and energy E are obtained in the 


usual way: 


dAo Of 
ee 
Fs dew oP OL Ja, 


fo 
B= At TS=N| fe-1(—) | 


(20) 


(21) 


oT 


The specific heats ¢., cp can also be expressed by 
differentiating E with respect to T either with 
a=const., or with p=const. using the equation 
of state (13). 

A calculation of the latent heat of melting* is 
not possible without a knowledge of the corre- 
sponding formulae for the liquid state which our 
theory does not provide. 

The thermodynamical properties of the crystal 
isolated or in equilibrium with its liquid can be 
expressed in terms of the functions fo, f1, fi, fie, faa 
of a and TJ and their derivatives. We proceed 
now to derive explicit expressions for these. 
For this purpose we have to make approximations 
of a different kind. 

The first approximation is the assumption that 
the thermic motion can be considered as harmonic. 
To justify this step we recall the simple model of 
Debye consisting of one particle moving in a 
nonlinear field of force. He uses this to explain 
the thermal expansion as the shift of the mean 
position due to the asymmetry of the amplitude 
in the nonharmonic field, and he shows that the 
vibration about this (shifted) mean position can 
be considered as harmonic. The situation be- 
comes much clearer in the case of a real lattice 
if we distinguish sharply between molar pa- 

*L. Brillouin, reference 8, says that a theory which 
considers melting as vanishing of shearing elasticity does 
not lead to a latent heat. This is quite correct insofar as 
an approach to any equilibrium from one side is insufficient 
for calculating the latent heat; but such a theory is not in 
contradiction with the existence of a latent heat (as, for 


instance, is the case for some types of so-called \-point 
transitions). 
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rameters—the dimensions of the cell (corre- 
sponding to the mean position in Debye’s model) 
—and the vibrating coordinates proper. Even if 
we assume the law of force for the latter to be 
linear, the averaged force acting on the molar 
parameters becomes nonlinear, and this suffices 
to explain thermic expansion. We shall show, 
moreover, that it also suffices to give a qualita- 
tive explanation of melting. Whether these 
formulae are accurate enough for quantitative 
purposes can only be judged after a thorough 
investigation which we shall carry out later. 
In any case, there is no difficulty in principle in 
developing the free energy for nonharmonic 
vibrations. 

For harmonic vibrations a well-known pro- 
cedure leads to the following expression for the 
free energy :8 


1 


| La y]| aa Le Pyei(e 0 — 1) 
> Gse'(e-*Ua) — 1) 
l 


where / is the index of the cell, representing 
really three integers /;, lo, 13; (la) =laitleas 
+lsa3, and ger’, gry’, °° the second derivatives 
of the potential energy between two lattice 
points in the relative position given by /(1;, /2, 13). 
The determinant (25) can be written as a 
triple sum; one finds after some calculations: 


[xy] =64>) & 


(I>0) (l’>0) (1’’>0) 


la) (Gia) _ (F''a) 
x {ll} eas ; sin? ; 








(26) 


where >- indicates that the sum is to be extended 
150 
only es half of the lattice points, in such a 
way that of two points J, J2, J; and —l1, —l2, —/s 
one is to be omitted. The coefficient {J/'l’’$ 
depends on the indices of 3 points, /(/iles), 
Vie ly), 0711s’), and can be written. in 
such a way that it is expressed by the invariants 
of the distance vectors r', ©”, r/’’ of these three 
points (where r’=/,a;+/2a2+/sa3), namely (r’)?, 


13See, for instance, M. Born, Atomtheorie des festen 


ipzi here 
Zustandes (Leipzig, 1923), p. 677. The method used 
is due to O. Sacre Ann. d. Physik (4), 51, 237 (1916). 


un 
oO 
wn 


hw 
A=%)+3NkT log —, (22) 
kT 
where po is the potential energy of the non- 
vibrating, but homogeneously deformed lattice, 


® a mean frequency (per 27 seconds) given by 
log @= g(log |[x, y]|)w—% log uw. (23) 


Here yu is the mass of the particles; the average 
is taken over the phases aj, a2, a3 of the waves, 


al f J flea, a2, as\derdasdas, (24) 


and |[x, y]| is the determinant of the secular 
equation for the waves with phases aj, ao, a3. 
For a simple lattice one has 





f= 


Lees" —1) Leal —1) Doeet(e-ilo 1) 
l l 


(25) 





(re rh) (r’-(r’ Xr)? We introduce the 


operator 


D=(1/r)(d/dr) (27) 


and write [Dg(r) ],-,.=Dg!'. Then we obtain: 
{W'1""} =D": Do” - De" +(t')?D?9'- De" Deo 
+4" X1r!"’)2. Deg!- D2g" +: Deg!” 


+4! (e!"X £!”))? D8g! Deg" 2g", (28) 
The invariants of the distance vectors depend 
only on the scalar products of the vectors 
aj, ae, a3, therefore, on account of (5), on the 
strain components x., -** ys, ***. Introducing 
(28) in (26) we get |[x, y]| as a function of the 
Lge ey ys; and of aj, a2, a3; averaging 
log |[x, y]| over these angles (as defined by 
(24)) leads to the value of @, (23), as a function 
of x;, «+: y3, °°*. This procedure, though quite 
definite and simple in principle, is in practice 
very involved and leads to formidable integrals. 
Therefore I shall use here a rather rough approxi- 
mation, which is suffcient to show that our 
theory represents the main features of lattice 
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If we assume g(r) as a power series in 1/r, 
beginning with 7” where the exponent m is 
large, we see that in the neighborhood of the 
minimum of ¢ the higher derivatives of ¢ will 
predominate. We proceed on the assumption 
that we can omit in (28) all terms except the 
last one which contains the highest derivatives 
and is proportional to the square of the volume 
of the parallelepiped determined by three lattice 
points: 


CU) 2 (rl. (r! Xrl'))?- Dg! D2g" - D2g"”’. (29) 


We assume further that the potential ¢(r) 
decreases so quickly with the distance that only 
the few next neighbors have to be taken into 
account. We shall—rather arbitrarily—take the 
first and second neighbors in the static terms 
(@o), but only the first neighbors in the vibra- 
tional terms (@). 

Finally, to fix our ideas, we consider one 
special case of a cubic lattice. The simple lattice 
is statically unstable if only the first neighbors 
are effective, and becomes stabilized only by 
the action of the second neighbors; this is the 
reason why this type practically never occurs in 
nature. The body-centered and the face-centered 
lattices occur in nature. I have decided (rather 
unfortunately, as was found later) to perform 
the calculations for the body-centered lattice since 
it has the coordination number 8 (instead of 12 
for the face-centered lattice), and this number 
of the first neighbors determines the number of 
terms in the sum (26). 

Then one obtains after some calculation 


I Tx, y || = 64 X16] a:acas|2- (A,Ri+A2Ro 


+4sRs+A.R,), (30) 


where A,, R, have the following meaning: The 
distances and phases for 4 of the 8 neighbors are: 


l=1(1 1 1) 11=d1+a2+a3 
20 1aiet) re= —d1+de+a3 
341-11) = r3=a1—a2+a3 
4111-1) = r4s=a;+a2—a3 (31) 


oF a)=aitartas, 

(2, aw) = —aitartasz, 

(3, a) =ai—az+az, 

(4, a) =a;tar—az. 
Then 
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(2, a) (3, a) | *(Apa) 
1=sin? sin? sin? ; 
2 
R,=D*¢2° D?¢3: D* ga, 
(3, a) (4,a) (1, a) 
Ag=sin? sin? sin? mer 
R2=D*¢3: D?¢4: D’ ¢1, (32) 
(4, a) (1, a) (2, @) 
Ag=<sill- sin? sin? , 
2 
R3=D?¢4: D?¢1- D* go, 
(1, a) Q;a)~ “Gia 
A,=sin? sin? sin? a. 


Ry=D? 91: D’ go D?¢3, 


where D?¢,=D?¢(r;) etc. 

The functions R, have to be developed with 
respect to the x,, --- yz, -::. For this purpose 
it is convenient to use the fact that any function 
of the distance, f(r), in which 7 is replaced by 
(r?+2p)*, has the following development : 


F?+-2p)*=f(r) E Dir) pe DY (1) =p ea 


just as a Taylor series where differentiation is 
replaced by the operator D. 

For the four neighbor points (31) we have 
11 =fo=r3=rs=ar/3 where a is the smallest 
coordinate difference in the lattice. The changes 
of r’ for a deformation x, +--+ yz, -+- are 


pr=A(xXetVyte-+y-+2r+%x,), 
p2=O7(xzt+yyt22+¥2—22—Xy), 
ps=Q(X2+Vy+2:—V:+22—Xy), 
ps=2(X2+Vy+3.—-V¥2-—S,+4x,). 


(34) 


With the help of (33) and (34) it is easy to 
develop the R,. The square of the determinant 
|aiacas! is equal to the determinant of the scalar 
products a,-a, and hence is expressible by the 
Xz, *** Vz, ***; One obtains for the terms up to 
the second order 


|aracas| =a3{1+(x,+y,+2.) —(y2+s2+x,2) 


+ (Yy2+22%2+%2yy) + es eee CS 


These developments must be substituted in (30) 
and then the averaging performed, as prescribed 
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in (23). We omit the very lengthy calculations and write down only the result :* 
to =4 9, +34 +3kT {> log (D293) —log T+loga+C}, 
D3 
f= —ya*p=4a*Do; +4a?Do, +er( 1440 =), 
D* 93 
Ds Diy5\* 
fa=yate, =4a'D¥es+160'D ge +kT| —2+ 404 =—( “) |} 
D? 93 D* 0; 
D4 D3 2 
fiz=ya%ey. =4a!D2 95 qe ta of ~) | 
D?o; \D*o3 
Dte; 161,.+5/7D3 
fss=a%ey, =4a'!D2¢, +47 —1+ 40 el =) |}: 
D? 93 9 D* 93 
where the constant C is given by between ;; and 4. We shall write 
25/3 [1/6 Iz=%—€, where 0<e<;%. (38a) 
C=log ae (36a) 
Ryu One can check the long calculations with the 


Here the lower indices, as ia gs, vu, Do; etc., 
indicate the square of the distance (in units a) 
for which the function is to be taken, e.g. 


gs=9(av/3), gsi=e(ar/4), 
1 deg(av/3) 


As dav/3 


According to our convention we have taken the 
first and second neighbors for the static terms, 
the first only for the kinetic terms. 

In these expressions some numerical integrals 
appear, J) and I». As Jy is only contained in the 
constant C its value is of no interest. The other 
constant is 


A, 2 
—f j || pear ot Sol i ee 
~ (Qn)3 (A, +A2+A3+A,)? (38) 


(37) 





I have not succeeded{ in finding an accurate 
value of J:; but it is easy to see that it lies 


*The formulae contain the well-known result that 
central forces, treated statically (T=0), lead to the 
Cauchy relation C12 =Ca4} but they show at the same time 
that for T>0 this is not the case. This is worth noting as 
deviations from the Cauchy relation are usually con- 
sidered as a criterion for noncentral forces which is not 
necessarily true. 

+ Note added in proof —T. du Cros has found by numer- 
ical calculation that Jz2=1.85/16, e=2.15/16, for the 
body-centered lattice. 


help of the formulae (13a) which are indeed 
satished (observe the meaning of Dos, Dy, -:: 
etc. as explained in (37)). 

The last step now consists in choosing a 
suitable function ¢(r). We take a law of the form 


nm LSTONG ASN 
(-2(7) CZ) 
n—-m MNT TONG Ts 
dg unm (C ny" (= “\") 
eae To N—-™M 


vanishes for r=7o; 7o is therefore the equilibrium 
distance for two particles. The frequency 10 of 
vibration of one particle about this equilibrium, 
if the other particle is fixed, is given by 


d*o u 
Aru? = -(=) =—mn. 
dr? r=ry 10" 


u= — o(fo) is the dissociation energy of a molecule 
formed of two particles, which can also be 
expressed as a temperature O0=u/k. We intro- 
duce further instead of 7p its projection dp on 
the cubic axes, given by 3ac’?=ro*, and the 
corresponding volume Vo=WNydo® which the 
crystal would have at zero temperature if only 
first neighbors interacted (the real volume at 
T=0 is a little smaller). We shall use further 





g=uUu 


then 
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a quantity of the dimensions of pressure, 
pbo=RO/Vo, where R=RN is the gas constant. 
Our law of force contains four constants, the 
numbers m, m and any two of the quantities 7 
(or do, or Vo), vo, #, O, Po, Which are connected by 


To? 47? P 
u=kO=—upvo"Po?, 


34/3. mn 
Po Vo — RO. 





Vo=Nyae=Ny 


(40) 


We shall use a final simplification by choosing 


n=2m. (41) 


Lennard-Jones! has found this assumption 
successful in the theory of the virial coefficient 
of gases, for the special value m=6, and he has 
used it also in his new order-disorder theory of 
melting, quoted above. 

If we divide Eqs. (36) by RO, the right-hand 
terms are dimensionless, independent of 79. On 
the left-hand side the coefficient of — p, ¢11, €12, C44 
becomes 


ya? ydo* a? (=) ~(- 3 
RO RO a)? RO Yo ney z=) , 

We choose now /9 as unit of length and replace 
the variable a=r/3 by 


1 1 m 
nbn Gay. 
Ses av/3 


Then the domain 1<r€ «, extending from the 
equilibrium position of the first neighbors to 
infinity, corresponds to the range 








(42) 





OS EX. (42a) 
As 
1 t 
—=——= (1-2) 
| Y ™ 
yg and all derivatives, Dg, gs, = * Do, <5 


are polynomials in £; for instance 


9 


a* r 


ne ei 2m 
5 amare thers 8 3: 


The result of a long series of reductions is the 


1s TE Lennard-Jones, Physica (IV) 10, 941 (1937). 
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following expressions : 
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TABLE I.* 








g V1—€}| Fil—) | Fo(é) |} Gié) | Ga(é) AMi(é) | He(€) | Hs(€) 























—0.25 || 1.118 || —5.00] 63.3 2.99 | —3.91 9.53 2.29] 1381 
—0.20 || 1.095 || —3.84] 56.3 3.00 | —5.54 9.73 3.26 | 136 
—0.15 |} 1.027 || —2.76] 49.7 3.00 | —7.03 9.96 4.47} 142 
—0.10 || 1.049 || —1.76] 43.4 2.99 | —8.37 10.24 6.00 | 150 
—0.05 |} 1.025 || —0.84] 37.5 2.96 | —9.56 10.58 8.00} 159 
0 1.000 0 32.0 2.93 | —10.60 11.00 10.67 | 171 
0.05 |} 0.975 0.76 | 26.9 2.88 | —11.49 11.53 14.37 | 186 
0.10 || 0.949 1.44) 22.1 2.83 | —12.23 12.22 19.72 | 207 
0.15 || 0.922 2.04 17.7 2.76 | —12.82 13.15 28.1 237 
0.20 || 0.894 2.56 | 13.7 2.68 | —13.26 14.50 42.2 285 
0.25 || 0.866 3.00} 10.0 2.59 | —13.56 16.60 68.9 367 
0.30 |} 0.837 3.36] 6.72 || 2.50 | —13.70 20.33 ; 539 
0.35 || 0.806 3.64] 3.81 || 2.39 | —13.69 28.8 341. 1054 
0.40 || 0.775 3.84} 1.28}) 2.27 | —13.54 67.0 | 2477. 5486 
0.45 || 0.742 3.96 | —0.88 || 2.14 | —13.23 || —73.0 | 4154. 6144 


*I have to thank Dr. E. W. Kellermann for computing these tables and for 
helping me in the subsequent numerical work. 


Table I contains the numerical values of these 
functions for m=6; this value corresponds to 
London’s quantum-mechanical explanation of 
van der Waals force and has been successfully 
used by Lennard-Jones, as mentioned above. 

The functions F;, F2 represent the statical 
effect of the first neighbors; Gi, Ge that of the 
second neighbors. As G; is not small compared 
with F,; the omission of all the more distant 
neighbors seems rather doubtful. Further, G2 is 
negative, therefore c1; is smaller than ¢ci2 and C44 
(at least for T=0, but really, as we shall see, 
throughout). Here it turns out that the choice of 
the body-centered lattice is unfortunate; for 
€11<¢12 means instability (of the type leading to 
a ‘‘gel’’ state). In fact, one sees without calcula- 
tion that for central forces acting only between 
first and second neighbors the face-centered 
lattice (densest packing) will be the most stable 
arrangement.* 

I have, nevertheless, calculated graphs for all 
physical quantities. (We may assume that the 
stability with respect to the condition ¢11>¢12 is 
maintained by other forces, not taken into 
account.) We shall concentrate our attention on 
C44 and investigate whether our thesis that ¢u:=0 
means melting gives reasonable results. 

It is easy to compute the quantities (43) as 
functions of p and T with the help of Table I, 
as the expressions (43) are linear in 7/0. I have 
used a graphical method in which these linear 
functions are represented by straight lines; it is 
obviously sufficient to calculate two points of 
each line from the table. From these systems of 
straight lines (one for each of the quantities 


*The calculation for the face-centered lattice will be 
performed by my collaborator, Mr. Teissier du Cros. 
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P, Cur, C12, C44) ONE Can construct curves which 
give &— (or the molar volume V) and the elastic 
constants as functions of p and T. 

Figure 1 represents the isobars of thermic 
expansion, £ (or V) for constant p as function of 
T/©. The V-scale is given in two forms; in the 
first one V is referred to Vo as unit, the quantity 
used in the theoretical considerations (repre- 
senting the volume of the lattice which it would 
have for T=0, p=0 if only the first neighbors 
interact); in the second scale V is referred to 
V(0), the actual (calculated) volume for T7=0, 
p=0. The figure shows that £ increases first 
linearly with 7/0, but later much quicker, until 
the increase becomes infinite. This extreme point 
is, of course, the absolute limit of stability; we 
have, nevertheless, continued the curves as 
dotted lines beyond these points. 

Figure 2 shows the effect of compression with 
the help of some isotherms é (or V) as function 
of p/po. 

Figures 3, 4, 5 and 6 represent the elasticity 
constants as functions of 7/@ for different 
hydrostatic pressures. As I could not determine 
the real value of the integral « I have calculated 
C44 for two values of e, namely e=0 (the lower 
limit) and «= 7s, which presumably is not far 
from the correct value. All constants decrease 
first linearly with 7/0, but turn downward at 
last and reach the zero line at steep angles; the 
values of 7/© where ¢44 vanishes should be the 
melting points for different pressures. We see 
that the behavior of the curves for ¢44 is rather 
different for e=0 and e=,;; in the latter case, 
Fig. 6, the zero points of c4, are not the maximum 
values reached by 7/0 on the curves. 

There is an arc on every curve on which to 
each value of ¢44/po there belong two values of 
T/@—a feature which could not well be foreseen 
from the general formulae (36). Those parts of 
the curves where C4, is not unique correspond 
apparently not to stable states of the lattices 
(or at least only to metastable, ‘‘overheated”’ 
states). The value of ¢4, jumps by a finite amount 
to zero at the melting point. Fig. 7 represents 
the melting curve p(T), as taken from the zero 
points of cys for e=0 and e=7¢; they do not 
differ much, but the latter is probably nearer to 
the truth. We see a slight increase of melting 
temperature with pressure. 
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If we now compare the curves of ¢11, C12, C44 
and particularly the zero points on them, we 
see that the curves of ¢i2 and C44, for «=O, 
practically coincide. The ci; curves lie lower and 
are flatter; this, as we have said already, is a 
consequence of the fact that the body-centered 
lattice is not the densest packing. 

We have now to consider whether the order of 
magnitude is right. This depends on the values 
of © and po. 

We introduce in (40) 


n Nyro? 
Vo= Nya? = 7 
cy 





M=,zN; 


then we get the dissociation energy u for a pair 
of particles in the form 


2 


ot Vo\! 
niat(—) : 
6 ¥ 16 N 


D=Nu=RO= 





(47) 


If we express this in kcal, we find with R=2 cal. 


D 
O=D-500°, py=—-41,400 atmos. (48) 
Vo 
For D=5, V)=20 eg. we have @©=2500°, 
y 
fae : “Yo) 
110 
7.09 re) 
I tea 
9.9Fbo 4 
£00 
°.9/ |-a2 
0,88] ~03 a9s 
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Po= 10,000 atmos. Our scale therefore covers the The expression (50) is Lindemann’s well- 
domain of high pressures investigated by known formula with a definite numerical con- 
Bridgman. stant. Lindemann’s empirical value is C=2.75; 

Figure 7 gives the melting point Tas some it is derived from those frequencies which allow 
numerical fraction of © for each value of p, say the representation of specific heat with the help 


Tne (49) of Einstein s formula. Our Vo cannot be directly 
considered as a vibration of the lattice. 
Introducing this in (47) and solving with respect One can calculate the Debye frequency 
to vo we get 
Vp =Cm(3.N/42V)3, (52) 


= 2 V3 2V7,—3 
pe tal) Ve where ¢, is a mean velocity of sound, from the 
= CX0.77 X10!2(7»,/MV>')?; (50) theoretical expressions of ¢11, C12, C14; but as the 
Debye formula of specific heat is valid only for 
deep temperatures, this calculation has nothing 
to do with our theory which is concerned with 


if we take for a the value from Fig. 7 which 
corresponds to p=0, a=0.35, we obtain 





Pree high temperatures (and pressures). I shall only 
Cat uit (51) show that yp differs from yo only by a numerical 
tT\ @ factor. As ¢11, C12, C44 are proportional to po we 
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have 

pCm? =Bpo, (53) 
where p is the density and 8 a numerical factor. 
Now p=M/V which we can replace by M/Vo 
(since, according to Fig. 1, V/Vo=0.93 for 
T=0). Then we have with the help of (47) 


Cn? = BVopo/M = BRO/M = (x°B/6 ¥ 16) ¥?(Vo/N)}, 


hence 
vp = 0.506? vo. (54) 
All numerical values calculated in this paper 
are valid only for lattices of the body-centered 
cubic type and for the assumed law of force with 
m=n/2=6; for substances of.this kind a “‘law 
of corresponding states’’ holds with respect to p 
and T if these are expressed as multiples of po 
and @. I shall illustrate this by a remark about 
the latent heat of melting which, as we have 
said, cannot be calculated from the theory, 
since the change of volume depends on the 
properties of the liquid state. But if AV is given 
we can apply the Clausius-Clapeyron formula 


Q=T(dp/dT)AV 
and take dp/dT from Fig. 7. We find 


dp/dT =11.1p9/0=11.1R/Vo, 

Q=11.1(AV/Vo)RT n= 32. (55) 
The latent heat is proportional to the melting 
temperature and to AV/V with a numerical 
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factor depending only on the lattice type and 
the exponent 7. 

I wish to make a final, quite hypothetical 
remark; it concerns the dotted parts of the ex- 
pansion curves of Fig. 1. They represent com- 
pletely unstable states of the lattice; but it can 
be expected that a more perfect theory which 
takes account of the disordered state of the 
liquid (as the theory of Lennard-Jones) will 
modify this part of the curves in the way indi- 
cated in Fig. 8. The ordinate of the melting 
point would intersect the curve in three points; 
the middle one would not represent a stable ° 
state, but the highest intersection point would 
correspond to the liquid in equilibrium with the 
solid. The liquid state would be represented by 
the continuation of the curve, exactly as in van 
der Waals theory of condensation. The ex- 
istence of a critical point between liquid and 
solid seems therefore not impossible. 

I expect that the method explained in this 
paper, properly improved, will be capable of 
accounting for the properties of strength, sliding 
and breaking of crystals. We have only to treat 
another stress component in the same way as 
we have here done with the pressure; we have 
to calculate the limit of stability with increasing 
temperature for a finite stress X, (breaking) or 
X,, (sliding). This kind of problem is of special 
interest in connection with Bridgman’s experi- 
mental investigations on breaking strength. 
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Bridgman has also discovered numerous allo- 
tropic transformations under high pressure; it 
will be interesting to see whether these can be 
accounted for by comparing the free energy of 
different lattices. 

Metals need a special treatment as the forces 
arising from the free electrons are not central. 

The extension of the theory to higher and 
deeper temperatures (anharmonic vibrations 
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giving terms with 7°, quantum effects terms with 
1/T) will also be considered. 

A complete account of the calculations will 
be published in collaboration with Dr. R. 
Schlapp. I have to thank him and Dr. K. 
Fuchs for valuable assistance in discussing the 
problems and checking the calculations, and Dr. 
E. W. Kellermann for computing the numerical 
table. 
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